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ABSTRACT 



After the approval by the Italian Space Agency of the LARES satellite, which 
should be launched at the end of 2009 with a VEGA rocket and whose claimed 
goal is a ~ 1% measurement of the general relativistic gravitomagnetic Lense- 
Thirring effect in the gravitational field of the spinning Earth, it is of the utmost 
importance to reliably assess the total realistic accuracy that can be reached by 
such a mission. The observable is a linear combination of the nodes of the existing 
LAGEOS and LAGEOS II satellites and of LARES able to cancel out the impact 
of the first two even zonal harmonic coefficients of the multipolar expansion of 
the classical part of the terrestrial gravitational potential representing a major 
source of systematic error. While LAGEOS and LAGEOS II fly at altitudes of 
about 6,000 km, LARES should be placed at an altitude of 1,450 km. Thus, it 
will be sensitive to much more even zonals than LAGEOS and LAGEOS II. Their 
corrupting impact Sfj, has been evaluated by using the standard Kaula's approach 
up to degree I = 70 along with the sigmas of the covariance matrices of eight dif- 
ferent global gravity solutions (EIGEN-GRACE02S, EIGEN-CG03C, GGM02S, 
GGM03S, JEM01-RL03B, ITG-Grace02s, ITG-Grace03, EGM2008) obtained by 
five institutions (GFZ, CSR, JPL, IGG, NGA) with different techniques from 
long data sets of the dedicated GRACE mission. It turns out 5fi ~ 100 — 1, 000% 
of the Lense-Thirring effect. An improvement of 2 — 3 orders of magnitude in the 
determination of the high degree even zonals would be required to constrain the 
bias to « 1 - 10%. 

Subject headings: Experimental tests of gravitational theories; Satellite orbits; 
Harmonics of the gravity potential field 
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Introduction 



In the weak-field and slow motion approximation, the Einstein field equations of 
general relativity resemble to the linear Maxwell equations of electromagntism. Thus, 
a gravitomagnetic field, induced by the off-diagonal components goi,i = 1,2,3 of the 
space-time metric tensor related to the m ass-energy currents of the source of the 



gravitational field, arises (IMashhoon 



particles and electromagnetic waves ( 



200 7j) . It affects in many wa ys the motion of test 



Ruggiero and Tartaglia 



2002|). Perhaps , the most 



famous gra vitomagnetic effects are the precession of the axis of a gyroscop e (IPugh 



Schifi 



I960) and the Lense-Thirringr 



precessions (ILense and Thirring 



1959 



19181 ) of the orbit of 



a test particle, both occurring in the field of a central slowly rotating mass like a planet. 
The measurement of the gyroscope precession in t he Earth's gravitat i onal fi eld has been 



the goal of the dedicated space-basedj GP-B mission (jEveritt et al 



197 



20011 ) launched in 



2004; its data analysis is still ongoing. The target accuracy was originally 1%, but it is now 
unclear if it will be possible to finally reach such a goal because of the onset of unexpected 
systematic errors during the mission. 

In this paper we critically discuss the possibility that the LARES mission, recently 
approved by the Italian Space Agency (ASI), will be able to measure the Lense-Thirring 
precession with an accuracy of the order of 1%. We recall that such a relativistic effect 
consists of a secular rate of the longitude of the ascending node Q of the orbit of a test 
particle 



c 2 a 3 (l 



l\3/2 ■ 



1 According to 



Pfisterl (120071 ). it would be more correct to speak about an Einstein- 



Thirring-Lense effect. 

2 See on the WEB http:/ /einstein. stanford.edu/ 
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where G is the Newtonian constant of gravitation, c is the speed of light in vacuum, S is 
the proper angular momentum of the central body, a and e are the semimajor axis and the 
eccentricity, respectively, of the test particle's orbit. 



In ( Van Patten and Everitt 



1976a 



van Patten and Everitt 



1976b! ) it was proposed 



to measure the Lense-Thirring precession of the nodes Q of a pair of counter-orbiting 
spacecraft to be launched in terrestrial polar orbits and endowed with drag-free apparatus 
to counteract the non-gravitational pertu rbations. A som ewhat equivalent, cheaper version 



of such an idea was put forth in 1986 by (jCiufolini 



19861 ) who proposed to launch a passive, 



geodetic satellite in an orbit identical to that of the existing LAGEOS satellite apart from 
the orbital planes which should have been displaced by 180 deg apart. LAGEOS was put 
into orbit in 1976, followed by its twin LAGEOS II in 1992; they are passive spacecraft 
entirely covered by retroreflectors which allow for their accurate tracking through laser 
pulses sent from Earth-based ground stations according to the Satellite Laser Ranging 
(SLR) technique. They orbit at altitudes of about 6,000 km (olageos = 12, 270 km, 
alageos ii = 12, 163 km) in nearly circular paths (e L AGEOs = 0.0045, e L AGEOs n = 0.014) 
inclined by 110 deg and 52.65 deg, respectively, to the Earth's equator. The Lense-Thirring 
effect for them amounts to about 30 milliarcs econds per year (mas yr" 1 ). The measurable 



Ciufolini 



(119861 ). the sum of the nodes of 



quantity was, in the case of the proposal by 
LAGEOS and of the new spacecraft, later named LAGEOS III, LARES, WEBER-SAT, 
in order to cancel to a high level of accuracy the corrupting effect of the multipoles of 
the Newtonian part of the terrestrial gravitational potential which represent the major 
source of systematic error (see Sectio n El). Although extensively studied by various groups 



(IRies et a] 



1989 



In ( llorio et al 



Ciufolini et al 



19981 ). such an idea was not implemented for many years. 



2002|) it was proposed to include also the data 



a different observable. Such an approach was proven in (llorio 



rom LA GEOS II by using 



2005al ) to be potentially 



useful in making the constraints on the orbital configuration of the new SLR satellite 
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less stringent than it was originally required in view of the recent improvements in 
our knowledge of the classical part of the terrestrial gravitational potential due to the 



dedicated CHAMP (http:/ /www. gfz-potsdam.de/pbl/op/champ/index_CHAMP.html) 



and, especially, GRACE (http://www.gfz-potsdam.de/pbl/op/grace/index_GRACE.html 



http: / /www.csr. utexas.edu/grace/ ) missions 



Reaching high altitudes and minimizing the unavoidable orbital injection errors is 
expensive; thus, it was explored the possibility of d iscarding; LAGEOS and LAGEOS II 



using a low-alt itude, near 



2006|), but in fllorio 



2002 



y pola r orbit for LARES (ILucchesi and Paolozzi 



2001 



Ciufolini 



2007al ) it was proven that such alternative approaches are not 



feasible. It was also su ggested that LARES would be able to probe alternati ve theories o f 



gravity (I Ciufolini 
2007bL 



20041 ). but also in this case it turned out to be impossible (jlorio 



2005bl . 



ASI recently made the following official announcement (http:/ /www. asi.it/SiteEN/MotorSearchFullTej 
"On February 8, the ASI board approved funding for the LARES mission, that will be 
launched with VEGAs maiden flight before the end of 2008. LARES is a passive satellite 
with laser mirrors, and will be used to measure the Lense-Thirring effect." The italian 
version of the announcement yields some more information specifying that LARES, designed 
in collaboration with National Institute of Nuclear Physics (INFN), is currently under con- 
struction by Carlo Gavazzi Space SpA; its Principal Investigator (PI) is I. Ciufolini and its 
scientific goal is to measure at a 1% level the Lense-Thirring effect in the gravitational field 
of the Earth. Concerning the orbital configuration of LARES, the ASI website says about 



VEGA that (http: //www.asi.it/SiteEN/ ContentSite.aspx?Area=Accesso+allo+spazio ): 
"[...] VEGA can place a 15.000 kg satellite on a low polar orbit, 700 km from the 
Earth. By lowering the orbit inclination it can launch heavier payloads, whereas 
diminishing the payload mass it can achieve higher orbits. [...]" In the latest com- 
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Ciufolinil (120081 ) writes that LARES 



munication to INFN, Rome, 30 January 2008, 
will be launched with a semimajor axis of approximately 7600 km and an inclination 
between 60 and 80 deg. More precise information can be retrieved in Section 5.1, pag 
9 of the document Educational Payload on the Vega Maiden Flight Call For Cube- 
Sat Proposals, European Space Agency, Issue 1 11 February 2008, downloadable at 



http://esamultimedia.esa.int/docs/LEX-EC/CubeSat%20CFP%20issue%201.pdf. It is 



written there that LARES will be launched into a circular orbit with altitude h = 1,200 
km, corresponding to a semimajor axis olares — 7, 578 km, and inclination % = 71 deg to 
the Earth's equator. Recent informationgj point towards a launch in November 2009 with 
an altitude oO h = 1,450 km, i.e. a = 7,828 km. Such an orbital configuration will make 



difficult to obtain a total accuracy of the order of 1%, as we will show in Section [21 



2. The major source of systematic error: the multipolar expansion of the 

Earth's gravitational potential 



In a realistic scenario the path of a probe is not only affected by the gravitomagentic 
field but also by a huge number of other competing orbital perturbations of gra yitational and 
non-g ravitational origin. The most important non- conservative accelerations (jMilani et al 



19871 ) are the direct solar radiation pressure, the Earth's albedo and various subtle thermal 



effects depending on the the physical 



state ( 



Inversi and Vespe 



Lucchesi et al 



2004 



1994 



Ries et al 



Vespe 



properties of th e sate l 



1999 



Lucchesi 



2001 



ite's surface and its rotational 



Lucchesi 



2002 



2003 



2004 : 



2003al ): however, the nodes of LAGEOS-like satellites are 



directly sensitive to them at a « 1% level only. Much more important is the impact that 



See on the WEB: http://www.esa.int/esapub/bulletin/bulletinl35/bull35f_bianchi.pdf 



4 D. Barbagallo, ESA-ESRIN, personal communication to the author, September 2008. 
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the oblateness of the Earth has on the satellite's dynamics. Indeed, the most insidious 
perturbations are those induced by the static part of the Newtonian component of the 
multipolar expansion in sphe rical harmon icaj Jg, I — 2, 4, 6, ... of the gravitational potential 



of the central rotating mass (IKaula 



19661 ): indeed, they affect the node with effects having 
the same signature of the relativistic signal of interest, i.e. linear trends which are orders of 
magnitude larger and cannot be removed from the time series of data without affecting the 
Lense-Thirring pattern itself as well. The only thing that can be done is to model such a 
corrupting effect as most accurately as possible and assessing, reliably and realistically, the 
impact of the residual mismodelling on the measurement of the frame-dragging effect. The 
secular precessions induced by the even zonals of the geopotential can be written as 



^geopot = £ ^ ^ 



(2) 



where the coefficients Q.e,£ = 2,4,6, ... depend on the parameters of the Earth (GM and 
the equatorial radius R) and on the semimajor axis a, the inclination i and the eccentricity 
e of the satellite. For example, for i = 2 we have 



3 (R 
■—n — 

2 a 



COS 2 



o2\2 



(3) 



n 



sjGMfo? is the Keple rian mean motion. They have been analytically computed up to 



20 in, e.g., (jlorio 



20031 ). Their mismodelling can be written as 



^opot < ^ I a 



1=2 



8Ji, 



(4) 



where bJi represents our uncertainty in the knowledge of the even zonal Ji of degree £; a 
way to assess them consists of taking the sigmas of the covariance matrix of a global Earth 
gravity solution provided that they have been realistically calibrated. It turns out that 



5 The relation among the even zonals Jg and the normalized gravity coefficients Ceo which 
are customarily determined in the Earth's gravity models, is Je = —\j2l+ 1 Ceo- 
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using the node of only one satellite is not possible because the gravitomagnetic signal would 
be swamped by the much larger mismodelled classical precessions. Thus, it is necessary 
to use other satellites to enhance the signal-to-noise ratio. The combination which will be 



used for me asuring the L ense-Thirring effect with LAGEOS, LAGEOS I 



obtained by 



Iorid ((2005a]) according to a strategy put forth by 



and LARES was 



Ciufolinil ril99fih ; it is 



SLAGEOS 



n 



SLAGEOS II 



+ C 2 Q 



S LARES 



, ci = 0.3586, c 2 = 0.0751. 



(5) 



and, by construction, is not affected by 5J 2 an d SJ4. The total Lense-Thirring effect, 
according to eq. ([1]) and eq. (j5J), amounts to 50.8 mas yr" 1 . The systematic percent error 
due to the mismodelling in the uncancelled even zonals J 6 , J 8 , ... can be conservatively 
evaluated as 



SfJ, < 



^LAGEOS _|_ Ci ^LAGEOS II _|_ ^^LARES 



50.8 mas yr 



-1 



100. 



(6) 



3. The impact of the high degree even zonal harmonics on a low-orbit LARES 



Nowadays, contrary to even a few years ago, a growing number of global Earth's 
gravity solutions, produced by different institutions by analyzing longer and longer data 
sets from GRACE with different approaches and techniques are becoming available. Thus, 
it is possible to obtain more robust and reliable evaluations of the total accuracy which 
could be realistically obtaine d with LARES acc ording to eq. (jSJ). Th e models we will 



use are EIGEN-GRACE02S dReigber et al 



from GFZ (Potsdam, Germany), GGM02S (ITaplev et al 



200 51) and E 



2007h from GSR (A ustin. Texas), ITG-Grace02s (IMaver-Giirr et al 



GEN- CG03C (iForste et al 



20051) and GGM 03S (ITaplev et al 



2005) 



( iMayer-Giirr 



20061 ) and ITG-Grace03 



20071) from IGG (Bo nn, Germany), JEM01-RL03B from JPL (NASA, USA) 



and EGM2008 flPavlis et al. 
USA). 



20081 ) from National Geospatial-Intelligence Agency (NGA, 
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Another very important point is that the low altitude of LARES forces us to consider 
much more even zonals with respect to the simpler original scenario for which a calculation 
up to degree I = 20 would have been more than adeq uate. Instea d, in all the p revious 
evaluations, the calculations were truncated at t = 20 (Iorio 



2005a ) or even less (jCiufolini 



20061 ). Thus, we numerically computed the node precession coefficients Qj of LAGEOS, 
LAGEOS II and LARES up to £ = 70 with two dif ferent softwares according to the 



Kaulal (119661 ). The figures determined in this 



standard and widely used approach by 
way have been successfully calibrated by compari ng them with those obtained from the 



coefficients analytically worked out by 



Ioriol (120031 ) up to i = 20; the same values have been 



obtained. As can be noted from an inspection of Table [TJ obtained by using the covariance 
sigmas for the aforementioned Earth's gravity models, the fact that LARES will be much 
more sensitive to the high degree even zonal harmonics than LAGEOS and LAGEOS II 
introduces a corrupting effect on the measurement of the Lense-Thirring precession which 
makes the expectation for a ~ 1% accuracy difficult to be implemented; it is true also in 
the cases in which only the formal, statistical sigmas were available. In Table Q] we quote 
both the upper bound obtained by linearly adding the absolute values of the individual 
mismodelled terms (SAV) and a more optimistic evaluation obtained by computing the 
square root of the sum of the squares of the individual biased terms (Root-Sum-Square, 
RSS); both the approaches yield very large figures. It turns out that improvements of at 
least 3 — 2 orders of magnitude in the even zonals, likely unattainable in the foreseeable 
future, would be needed to constrain S/j, at the ~ 1 — 10% level. It would be desirable that 
other researchers will independently repeat the calculations with different computational 
approaches. 
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Table 1: Systematic percent error Sfj, (%) according to eq. (jSJ) up to £ = 70. Second column: 
linear sum of the absolute values of the individual biased terms (SAV). Third column: square 
root of the sum of the squares of the individual biased terms (Root-Sum-Square, RSS). The 
orbital parameters of LARES are a = 7, 828 km, % = 71 deg, e = 0. For 5Jg, i = 2, 4, ...70 the 
sigmas of the covariance matrices of EIGEN-GRACE02S (calibrated), EIGEN-CG03C (cal- 
ibrated), GGM02S (formal), GGM03S (calibrated), ITG-Grace02s (formal), ITG-Grace03 
(formal), JEM01-RL03B (formal) and EGM2008 (calibrated) have been used in both cases. 



Model 5(i (up to I = 70: SAV) tf/x (up to I = 70: RSS) 



EIGEN-GRACE02S 


4, 443% 


2, 563% 


EIGEN-CG03C 


2, 177% 


1,256% 


GGM02S 


3, 595% 


2, 062% 


GGM03S 


911% 


524% 


ITG-Grace02s 


63% 


36% 


ITG-Grace03 


58% 


33% 


JEM01-RL03B 


235% 


135% 


EGM2008 


1,307% 


754% 
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4. Conclusions 

In view of the recent approval of the LARES mission, to be launched at the end of 2009 
with a VEGA rocket into a circular orbit with semimajor axis a = 7, 828 km and inclination 
i = 71 deg, we discussed the possibility that the LAGEOS, LAGEOS II and LARES 
satellites will be able to measure the Lense-Thirring effect at an accuracy of the order of 
1%. The observable will be a linear combination of the nodes of the three spacecraft which, 
in fact, should be affected by the non-conservative perturbations at a « 1% level. A major 
source of systematic bias is represented by the uncertainty in the static part of the even 
zonal harmonic coefficients of the multipolar expansion of the Newtonian component of the 
terrestrial gravitational potential. The low altitude of LARES, 1,450 km with respect to 
about 6,000 km of LAGEOS and LAGEOS II, will make its node sensitive to much more 
even zonals than its two already orbiting twins; it turns out that, by using the sigmas of 
the covariance matrices of some of the latest global Earth's gravity solutions based on long 
data sets of the dedicated GRACE mission, the systematic bias due to the mismodelled 
even zonal harmonics up to £ = 70 will amount to ~ 100 — 1,000%. An improvement of 
3 — 2 orders of magnitude in the determination of the high degree even zonals would be 
required to reduce the systematic error to ~ 1 — 10% level. 
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